Electron correlations and bond-length fluctuations in copper oxides: from Zhang-Rice 

singlets to correlation bags 
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We perform first principles, multiconfiguration calculations on clusters including several CuOe 
octahedra and study the ground-state electron distribution and electron-lattice couplings when 
holes are added to the undoped d 9 p 6 configuration. We find that the so-called Zhang-Rice state 
on a single CuCU plaquette is nearly degenerate with a state whose leading configuration is of the 
form Cu d 9 - Op 5 - Cu d 9 . A strong coupling between the electronic and nuclear motion gives rise to 
large inter-site charge transfer effects for half-breathing displacements of the oxygen ions. Under 
the assumption of charge segregation into alternating hole-free and hole-rich stripes of Goodenough 
QLQi our results seem to support the vibronic mechanism and the traveling charge-density wave 
model from Refs.0, for the superconductivity in copper oxides. 



I. INTRODUCTION 

It is widely believed that the lattice degrees of freedom 
should play some role in the high-T c superconductivity 
in doped copper oxides. In a recent series of articles, see 
for example 0, 0, Q, J. B. Goodenough and coworkers 
predict for these systems the existence of strong lattice 
instabilities associated with two distinguishable (equilib- 
rium) Cu-0 bond lengths and suggest that long range 
"ordering" of bond-length fluctuations would stabilize 
some kind of traveling charge-density wave (CDW) in 
the superconducting state. Goodenough et al. propose a 
model where the hole-hole Coulomb repulsion and rela- 
tively strong antiferromagnetic (AFM) d 9 -d 9 couplings 
give rise at x — 1/6 doping and low temperatures to a 
configuration with alternating, 1/3-doped and hole-free, 
Cu-O-... stripes. Within the 1/3-doped Cu-O-... rows 
the holes occupy Cu-O-Cu bonding orbitals. The super- 
conductivity is associated in Refs. [HQ. El with itinerant 
heavy vibrons that originate from the coupling between 
some particular (dynamic) ordering of these Cu-O-Cu 
"hole bags" and in-plane phonons. 

The present paper is an attempt towards a better un- 
derstanding of the nature of the local electron-electron 
and electron-lattice interactions in copper oxides. We 
employ ab initio methods from traditional quantum 
chemistry. Our calculations are performed on finite clus- 
ters including several CuOg octahedra. We pay special 
attention to the situation where one hole is added to 
the undoped, formally Cu d 9 Op 6 , configuration. Results 
from a square 5-octahedra cluster indicate that the so- 
called Zhang-Rice (ZR) singlet on a CUO4 plaquette 3 
may become "unstable" with respect to a lower symme- 
try oxygen- hole state with a Cu d 9 - O p 5 - Cu d 9 leading 
configuration. This is just the 1-hole Cu-O-Cu correla- 
tion bag of Goodenough jlj . We further investigate local 
interactions within a 1-hole, 4-octahedra linear cluster. 
We find evidence for strongly anharmonic, Jahn- Teller 
(JT) type effects for certain Cu-O bond distortions. The 



oxygen-atom vibrations are coupled with strong inter-site 
charge transfer (CT) effects. Our findings seem to sup- 
port a mechanism of the type proposed by Goodenough 
et al. P, 0, Q for the high temperature superconductiv- 
ity. 

It was suggested in Ref. Q that at x = 1/6 doping 1- 
hole Cu-O-Cu units could cluster to form spin-paired 
2-hole bags including four Cu centres. A hypothetical 
ordered configuration of such 2-hole bags is shown in 
Fig. 1(a). In a Zhang-Rice picture, the 1/3-doped stripes 
can be represented as in Fig. 1(b), for example. Our anal- 
ysis is carried out on 1-hole clusters corresponding to a 
somewhat intermediate arrangement, as we shall discuss 
below. 
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FIG. 1: (a) 2-(doped-)hole, 4-Cu correlation bags in the CuO 
plane at 1/6-doping as suggested by Goodenough |2J, see text. 
Only the Cu sites are shown. A hole-free (d 9 p 6 ) chain is 
drawn with open circles, (b) Hole-free and 1/3-doped stripes 
in a Zhang-Rice picture. The squares represent ZR singlets 
on Cu04 plaquettes. A relative shift between neighboring 
doped stripes should minimize the Coulomb repulsion. Intu- 
itively, a 2-hole correlation bag forms by "clustering" of two 
ZR singlets. 
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II. LOCAL CORRELATIONS IN THE CUO 
PLANE: CUO4 ZHANG-RICE SINGLETS 
VERSUS CUd 9 -Op 5 -CUd 9 CONFIGURATIONS 

First, we perform electronic structure calculations on a 
square cluster including a central CuC>6 octahedron and 
the four adjacent in-plane octahedra. This [CU5O26] clus- 
ter is embedded in a matrix of point charges that repro- 
duce the Madelung potential associated to the undoped 
La2CuC>4 crystal. The nearest neighbor cations are rep- 
resented by effective ion potentials. We use the tetrago- 
nal crystal structure of Lai.gsSrg .15CUO4 as reported by 
Cava et al. p|, with in-plane Cu-0 distances of 1.89 A. 

We rely on an ab initio, wave-function-based multi- 
configuration (MC) approach G]. Our starting point is 
a Hartree-Fock self-consistent-field (SCF) calculation for 
a closed-shell Cud 10 Op 6 configuration of the [Cu 5 2 6] 
cluster. We remove then six electrons. We construct 
thus a multi-determinant (or multi-configuration) wave- 
function where "on top" of a number of closed electron 
shells six electrons/holes are distributed in all the possi- 
ble ways over six orbitals. Such a configurational space is 
referred to as a complete active space (CAS) @. Those 
electrons and those orbitals defining the MC space are 
called active. Due to the relatively large size of the clus- 
ter we use effective core potentials (ECPs) for the Cu 
Is, 2s, 2p, 3s and the four bridging-oxygen (0£, see Fig.2) 
Is core shells. The ECP basis sets of Huzinaga et al. Q 
were applied for these ions, with the following contrac- 
tions: Cu 9sQp6d/2s2p3d and O 5s6pld/2s3pld. For the 
other oxygens we applied the basis set of Widmark et al. 
P| contracted to 2s2p. All calculations were performed 
with the MOLCAS program package . 

Unexpectedly, for an undistorted structure with identi- 
cal Cu-0 distances, the CAS SCF ground-state optimiza- 
tion converged to a broken-symmetry triplet state where 
the largest weight corresponds to a Cu c d^ 2 _ y2 - Oj p y - 
Cul q d x 2 _ y 2 configuration (see Fig. 2(a)), although our 5- 
octahedra cluster has full D^h symmetry ^(j . The 0\ p y 
and Cu* d x i_ y i atomic orbitals form in this geometry 
bonding and antibonding combinations. Contributions 
from the Cu c d x 2_ y 2 and p x /Py orbitals on ligands adja- 
cent to the Cu^ cation to this pair of bonding and an- 
tibonding orbitals, while also present, are much smaller. 
Due to correlation, not only the bonding orbital is occu- 
pied but a partial charge transfer into the antibonding 
orbital occurs. The latter has an occupation number of 
0.27, in comparison to 1.73 for the former. The forma- 
tion of such a C\i sq d x 2_ y 2-0\py bond favors antiferro- 
magnetic d x 2_ y 2~d x 2_ y 2 interactions over the other three 
bridging oxygens. 

It turns out that the so-called ZR state, where the oxy- 
gen hole is equally distributed over four anions (Fig. 2(b)), 
is approximately 80 meV above the broken-symmetry 
state with overall triplet spin multiplicity. Our data 
also indicate that it costs 44 meV to couple ferromag- 
netically those three Cu sg d x 2_ y 2 electrons to the Cu c 
d x 2_ y 2 (Fig. 2(c)). Next, we determine the positions of 




FIG. 2: (a) Leading configuration for the broken-symmetry, 
lowest-energy state in the geometry with identical Cu-O dis- 
tances, see text. A CUO4 plaquette and four in-plane Cu 
neighbors are shown, (b) Schematic representation of the ZR 
state on our cluster, (c) Leading configuration for the lowest 
excited state in the undistorted structure. The d-d couplings 
are ferromagnetic. 



the bridging oxygens 0\ that correspond to a minimum 
of the cluster total energy. According to our results this 
is the structure where the four Cu c -O h bonds are all 
about 5% shorter. For this arrangement, the ground 
state is a ZR like state with full D±h symmetry. No 
broken-symmetry solution could be obtained for this ge- 
ometry. Now, if we compare Fig. 2(b) with Fig. 2(c) (or 
Fig. 2(a)) we can imagine the following scenario: starting 
from the minimum-energy structural configuration with 
shorter Cu c -O h bonds and a ZR electronic ground state, 
(half-) breathing displacements of the Ob oxygens may 
result in a broken-symmetry state where the leading con- 
figuration is of the form Cud X 2_ y 2-Op y -Cud X 2_ y 2 and 
the oxygen 2p hole has the largest weight (i.e. is partially 
transferred, about 0.3 of an electron charge) onto a single 
anion. It is important to mention that the wave functions 
and the relative energies for the states discussed above 
remain largely unchanged when the active space is ex- 
tended to include more ligand 2p orbitals, such as the 
other linear combinations involving Ob p y /p x bridging 
orbitals. 

It would be useful to determine how the cluster charge 
distribution changes for oxygen-atom distortions on an 
adjacent CUO4 plaquette. However, due to hardware lim- 
itations we could not employ sufficiently accurate basis 
sets for those oxygens beyond the nearest neighbors of 
the central Cu ion. To study such effects, the "cheap- 
est" option is a 1-hole, 2-octahedra cluster. Calculations 
on this cluster reveal strong charge fluctuations between 
the two symmetric, minimum-energy geometrical config- 
urations where the in-plane Cu-0 distances on only one 
of the plaquettes are 5 — 6% shorter. For each of these 
distorted configurations, the doped hole has the largest 
weight on a single O4 plaquette and forms a ZR type sin- 



glet with the Cu 3d hole. In the undistorted geometry, 
the oxygen hole is mostly localized on the ligand bridg- 
ing the two octahedra. The energy of the high-symmetry 
state is about 200 meV higher as compared to the ZR like 
configurations. The results indicate strong anharmonic 
effects, i.e. a one-dimensional double- well, when shifting 
the bridging oxygen along the y axis. This in analyzed 
in the next section. 

For an isolated Cu c 06 octahedron (or an isolated 
5-octahedra square cluster) with no distortions, the 
simplest way to view the I-doped-hole ground-state 
wave-function is as a superposition of the four broken- 
symmetry, degenerate Cu^+-0~ solutions, see for ex- 
ample the discussion in Nevertheless, the calcula- 
tions on the 2-octahedra cluster show that there is a finite 
probability for the ZR singlet to hop to an adjacent pla- 
quette via such an intermediate broken-symmetry state 
with a dominant d?-p 5 -d 9 configuration. (Meta)stable 
d 9 -p 5 -d 9 configurations could actually be induced by 
longer range Coulomb interactions and AFM spin corre- 
lations in the doped plane. It seems then that the ex- 
istence of l-(doped-)hole Cu-O-Cu entities that could 
eventually cluster to form multi-hole correlation bags, as 
suggested by Goodenough is not an unrealistic hy- 
pothesis. 
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FIG. 3: (a) Possible charge-ordered configuration with 1- 
(doped-)hole, 3-Cu spin-singlet "correlation bags" along the 
1/3-doped stripes, see text. Only the Cu sites are shown. 
A hole-free chain is drawn with open circles, (b) Four con- 
secutive CUO4 plaquettes along the 1/3-doped chain. The 
CAS SCF calculations were carried out on such clusters. The 
O 2p hole is represented with black or gray circles. For cer- 
tain phases of the O-atom half-breathing vibrations the hole 
can hop along the chain. From left to right, the Cu and the 
bridging O atoms along the CuO chain are labeled as Cuu, 
O l b , Cu;, Ob, Cu r , Oft, and Cu rr . 



III. 1-HOLE, 4-OCTAHEDRA LINEAR 
CLUSTERS: BOND-LENGTH FLUCTUATIONS 
AND INTER-SITE CHARGE TRANSFER 

Inspired by the work of Goodenough et al. 0, 0, Q 
we investigate local many-body effects and electron- 
lattice interactions within a l-(doped-)hole, 4-octahcdra 
[CU4O21] linear cluster. We apply Ci v symmetry re- 
strictions, with xy and zy mirror planes (see Fig. 3), 
and all-electron basis sets: Cu 21sl5pl0d/5s4p3ci [l2fl. 
O 14s9p4<i / 4s3pld for the bridging 0& oxygens and O 
14s9p/4s3p for the other oxygens y|. 

Starting from an undistorted structure with identical 
Cu-0 distances, we perform a number of test calcula- 
tions to identify the geometry that minimizes the cluster 
total energy. However, only the chain Cu and Oh ions are 
allowed to move, along the y axis. We do not modify the 
positions of the in-plane oxygens adjacent to the CuO 
chain, O l a , or the positions of the apical ligands. We em- 
ploy a 5 electron / 5 orbital CAS, with five electrons being 
removed from the closed-shell Cu d 10 O p 6 configuration. 
It is found that the lowest-energy geometry corresponds 
to a distorted configuration like that shown in Fig. 3(b). 
The atom is shifted to the left by 6% of the high- 
symmetry Cu-0 bond length and the O^ atom is shifted 
by 6% to the right. For this geometry an oxygen 2p 
hole/electron is bound to the Cu; d x 2_ y 2 hole/electron 
into a ZR like singlet state. We give in Table I the com- 
position of the bonding (B) and antibonding (AB) active 
(natural) orbitals 01 on the Cu;04 plaquette, their oc- 
cupation numbers and the Mulliken populations (MPs) 



for the relevant O 2p and Cu 3c? atomic orbitals. The 
other three active orbitals have essentially Cun , Cu r and 
Cu rr d x 2_ y 2 character and occupation numbers of nearly 
1 . The overall spin multiplicity is doublet, with the Cu r 
and Cu rr 3d electrons coupled antiferromagnetically. 

Now, from this geometrical configuration we shift the 
O^ and Og nuclei back to the high-symmetry, undistorted 
arrangement in steps of 1% of the Cu-0 bond length. 
The orbitals from the previous step are used each time 
as input orbitals 0|. For the undistorted geometry the 
calculation converges to a state where the doped hole 
is mainly localized on the bridging O l b oxygen. Details 
regarding the composition and the occupation of the rel- 
evant orbitals are given in Table I. It was found that the 
energy of this O^-hole state in the undistorted geometry 
is 320 meV higher then for the ZR like state in the con- 
figuration where the Cu-O^ and Cu-Og bonds are 6% 
shorter. In agreement with the results from the previ- 
ous section, between these two configurations a transfer 
of about one third of an electron/hole occurs from the 
O l b p y to the 0„ p x , 0£ p y , and Cu; d x 2_ y 2 orbitals, see 
Table I. The total cluster energy is plotted for different 
positions of the O^ and 0£ oxygens in the left part of 
Fig.4. 

For systems with competing valence configurations like 
our hole-doped cluster, different sets of input orbitals 
can lead to different solutions. We performed calcula- 
tions where we forced the localization of the doped hole 
on the central 0£ p y orbital, by adding some extra pos- 
itive charge within our point charge embedding. Using 
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TABLE I: Contributions of various oxygen 2p and metal 3d atomic orbitals to the relevant Cu-0 bonding (B) and antibonding 
(AB) active (natural) orbitals [l3T |. In each case, the other three active orbitals have essentially Cu d x 2_ y2 character and 
occupations of nearly 1 . Coefficients smaller than 0.10 are not shown. Natural-orbital occupation numbers and Mulliken 
populations (MPs) are also given. Of, are bridging oxygens along the CuO chain, the O a ions form parallel, adjacent in-plain 
O chains, see Fig. 3(b). 
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a The fact that the doped holes have the largest weight on O 
Px/Py orbitals pointing to the Cu sites wa s also found in previous 
ab initio wave-function-based calculations ll5lll6l[T^ . Tendencies 
to ward s localization and broken-symmetry solutions were discussed 
in 0. 

b The chain Cu;— O;, bonds are 6% shorter. 



this output as a new set of starting orbitals for CAS SCF 
calculations with the initial embedding and no distor- 
tions, we were able to obtain a solution with a dominant 
contribution from the ...-Cu; d l 2 2 - Of p}- Cu r d\ 2 - 

L x — y o^y x —y 

... configuration. The energy of this state is 350 meV 
above the minimum corresponding to the ZR state on the 
distorted CU/O4 plaquette (see Fig. 4). For the O b -hole 
state, the Cu r d 1 2 2 -Cu rr d 1 2 2 AFM interactions are 

1 ' x — y 11 x — y 

not perturbed by the presence of the oxygen hole and 
its relative energy is somewhat lower, 320 meV. The fact 
that the solution with the doped hole mainly localized on 
the 01 ligand is less accessible can be understood on the 
basis of this energy difference between the two Ob-hole 
states. Regarding the data in Table I and Fig. 4 we also 
note that the Og-hole doublet wave-function calculated 
for the undistorted geometry misses full Di n symmetry. 
Attempts to converge to a fully symmetric wave-function 
were unsuccessful. 

The structural configuration with 6% shorter chain 
Cuj-0 bonds is symmetry equivalent to the situation 
where the chain Cu r -0 bonds are 6% shorter. These 
two distorted configurations represent the minima of a 
one-dimensional double-well potential. From one mini- 
mum to the other, the sum of the Mulliken populations 
of the Ob p y , O a Px , and Cu d x i_ y i orbitals on each of 
the Cu;04 and Cu r 04 plaquettes changes by about 0.7 e. 
The total Mulliken electronic charge on each of these pla- 
quettes changes between the two minima by about 0.5 e 
(not shown in Table I) 0] . It is well known that the MPs 
are basis set dependent. Still, variations in the charge 
population of a whole CUO4 plaquette (or CuO@ octahe- 



dron) is a reliable quantitative indicator. An illustrative 
result is also the variation of the total Mulliken electronic 
charge for the Cu/Oq and G\\uO§ octahedra when going 
from the state with a Cu^ d^ 2 _ y 2~O l b p y - Cu/ d\ 2 , lead- 
ing configuration to the valence structure where the hole 
is mainly on the 0\ oxygen. This value is 0.95 e. About 
the same amount of charge, one hole/electron, should 
be transferred along the chain between "traveling" next 
nearest neighbor ZR singlets. 

As mentioned in the previous section, we performed 
extra calculations on a 1-hole, 2-octahedra cluster. The 
doped hole is confined to a smaller region in space in this 
situation. In addition, charge relaxation effects within 
the neighboring chain octahedra can not be accounted 
for. Therefore, the CUO4-CUO4 CT is weaker when dis- 
torting the CuO bonds, about 0.3 e. For this smaller 
cluster it is possible however to use a more flexible active 
space, including all O 2p and Cu 3d orbitals. This cor- 
responds to distributing 3 holes over 43 orbitals. With 
such an active space we found a CUO4-CUO4 CT of 0.4 e 
for the half-breathing distortions. 

The results reported here are quite remarkable. They 
show that within a 1-doped-hole cluster including several 
CuOg octahedra, complex electron correlations give rise 
to quasi-localized O 2p hole states and strong charge in- 
homogeneity. Several competing O-hole valence config- 
urations were identified. The near-degeneracy between 
different hole configurations and the electron-lattice in- 
teractions cause strongly anharmonic, pseudo JT effects 
and facilitate the hopping of the charge carriers. Starting 
from a minimum-energy geometry where shorter "chain" 
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FIG. 4: Quasi-localized O-hole states on a 1-doped-hole, 4- 
octahedra linear cluster, see text. The energy maxima corre- 
spond to states with dominant Ot-hole character in the undis- 
torted structure. The two minima are related to ZR like states 
where the Cu-Ot distances on the Cu;04 plaquette (the left 
hand minimum) or the Cu-Ot distances on the CurCU pla- 
quette (the right hand minimum) are 6% shorter. From each 
minimum-energy structural configuration, the two Ob ligands 
are shifted back to the undistorted geometry in steps of 1% 
of the high-symmetry Cu-0 bond length. The energies of 
these states are shown with black dots. The lines are a guide 
for the eye. The last point on each of these curves corre- 
sponds to a configuration where the Cu-Ot, distances on the 
Cui;04/Cu r7 .04 plaquette are shorter, by 1%. Other calcu- 
lations were carried out for the case where the Cu-Of, bonds 
on the CuiOm and Cu r 04 plaquettes are gradually shortened 
starting from the undistorted structure and using the orbitals 
of the Oj-hole state as starting orbitals. Due to convergence 
problems, only few points could be obtained. These points 
are shown as black triangles. 



Cu-0 bonds on a Cu04 plaquette give rise to a quasi- 
localized ZR like state, a half-breathing motion of the 
chain oxygens with a phase difference of n between ad- 
jacent plaquettes (see Fig. 3(b)) leads to a CUO4-CUO4 
charge transfer of about 0.5 e over an energy barrier of 
350 meV (175 meV per chain oxygen). The minima of the 
one-dimensional double-well imply in-line O-atom dis- 
placements about the middle positions with amplitudes 
of only 0.11 A (±6% of the Cu-0 bond length). Although 
we did not investigate the effect of further Cu-Of, dis- 
tortions, it seems that a whole electron charge moves 
along the chain when a ...-CUO4-... ZR like singlet 
is "transferred" to a next nearest neighbor plaquette. 

It is possible that at about 1 / 6-doping an ordered ar- 
rangement of hole- free and 1/3-doped stripes is formed 
within the CuO plane pol | . Under this assumption our 1- 
hole, 4-octahedra cluster is related to the unit cell of the 
charge-ordered 1/3-doped stripe shown in Fig. 3(a) (such 
a "unit cell" actually contains only three in-line octa- 
hedra). We speculate then that longer range hole-hole 
Coulomb interactions associated with cooperative, corre- 
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FIG. 5: Possible representation of the traveling charge- 
density /spin-density wave along a 1/3-doped CUO4-CUO4- 
... chain. The O 2p holes are shown with black or gray 
circles. For each sequence, the spins on the first four pla- 
quettes display schematically fl9| the CAS SCF results. The 
ellipses indicate spin-singlet couplings. Dashed arrows indi- 
cate a "spin-flip" process. 



lated oxygen-atom displacements like those described in 
the previous paragraph along a chain of several Cuff- O^- 
Cui-Og-Cu r -0£ units should lower the energy barrier 
corresponding to the inter-site hole hopping and make 
possible an ordered, collective motion of the doped holes. 
Under the assumption of an alternating 1/3-doped/hole- 
free stripe configuration, our results seem to support thus 
a vibronic mechanism like proposed by Goodenough 0,0 
for the superconductivity in cuprates. It still needs to 
be investigated whether (and why) nearest neighbor 1- 
hole correlation bags would "cluster" along the chain to 
form a (traveling) CDW with spin-paired 2-hole, 4-Cu 
units separated by two antiferromagnetically-coupled Cu 
d? ions as suggested in Ref. Q (see Fig. 1(a)). The longer 
range interactions responsible for the formation of such 
paired oxygen holes should be however weaker than those 
causing the local JT/CT effects we illustrate in Fig. 4 and 
Table I. They probably merely modify the amplitudes of 
the O-atom displacements that determine the double- well 
potential. It should be pointed out at this stage that the 
CAS SCF results suggest that nearest neighbor oxygen 
holes have already antiparallel spins, see Fig. 5 . 

Obviously, reliable estimates for the strength and the 
effects of the longer range couplings, intra- and inter- 
chain, are required for quantitative predictions. We note 
that in this model the effect of the inter-chain Coulomb 
interactions should be stronger for CuO multilayer struc- 
tures like in the mercury compounds and might con- 
tribute to a higher T c . The critical temperature T c 
should be related to thermal fluctuations that disrupt 
the phase coherence among the kind of O-atom displace- 
ments described above. At temperatures well above T c 
the striped arrangement is probably destroyed. For a cer- 
tain temperature interval above T c , the transition from 
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a charge-inhomogeneous phase with parallel chains (or 
fragments of chains) of 1-holc, 3-octahedra units but 
"uncorrelated" bond-length fluctuations to a totally dis- 
ordered charge configuration could provide an explana- 
tion for the (71", 7r) — (7T, 0) transfer of spectral weight in 
the angle-resolved photoemission spectroscopy (ARPES) 
data I23. 

Our [CU4O21] linear cluster seems to be a reason- 
able choice for studying local correlations for a 1/3- 
doped chain of octahcdra. The quasi-localized char- 
acter of the oxygen hole (see Table I) justifies such a 
material model. This cluster model would probably 
be inappropriate for the 1 /2-doped chain (or stripe) in 
Lai.48Ndo.4oSro.i2Cu04 j24j, where the doped-hole or- 
bitals have much larger overlap. Also, although the 
CAS SCF method can not yield highly accurate ener- 
gies, such an approach is sufficient for an accurate de- 
scription of the electron distribution provided the proper 
active space is chosen. The near-degeneracy and CT ef- 
fects illustrated by our CAS SCF results should be then 
a genuine characteristic of the CUO4-CUO4-... chain. 
On the other hand, the CAS SCF data normally gives a 
rather qualitative (or semi-quantitative) description for 
potential-energy landscapes like that associated with the 
O-atom half-breathing displacements. A more elaborate 
treatment, e.g. multi-reference configuration-interaction 
or multiconfiguration perturbation theory is needed 
for calculating accurate potential curves and surfaces. 
Such calculations will be documented in an upcoming 
publication. 

A matter of concern with our calculations may be the 
fact that we used experimental structural data. It would 
be desirable to start, for example, from a fully opti- 
mized structural configuration of the undoped system. 
For such an optimized geometry, the in-plane Cu-Cu dis- 
tance would be of particular interest. It is known that the 
ab initio periodic Hartree-Fock (HF) calculations in the 
linear-combination of atomic orbitals (LCAO) formula- 
tion slightly overestimate the lattice parameters in tran- 
sition metal oxides, for example by 2% in MnO and NiO 
|2H and almost 4% in Cu 2 Hj. To reduce these devia- 
tions, post-HF techniques are needed [23. Nevertheless, 
structure optimizations for complex periodic systems like 
La2Cu04 at a correlated, post-HF level are not possible 
yet. Intuitively, we expect that for slightly larger Cu-Cu 
distances in our cluster the doped hole is somewhat more 
localized on the distorted CUO4 plaquette (see Table I). 
In this case the half-breathing O-atom vibrations should 
imply more CT between adjacent plaquettes but over a 
higher energy barrier. 

IV. CONCLUSIONS 

We study local electron correlation effects for hole- 
doped CuO clusters including several Cu06 octahedra. 
We use the structural data reported for the tetragonal 
lattice of Lai.85Sr .i5CuO4. Results of ab initio MCSCF 



calculations on a 1-hole, 4-octahedra linear cluster indi- 
cate a double-well potential for half-breathing O-atom 
displacements along the O-Cu-O-Cu-0 row. The Cu- 
O bond-length fluctuations are coupled with large inter- 
site charge transfer effects. For distortions of ±6% of 
the Cu-0 bond length, charge fluctuations of about 0.5 e 
are observed between nearest neighbor CUO4 plaquettes. 
It seems also that a whole electron charge would move 
along the chain when a ...-CUO4-... ZR like singlet is 
transferred to a next nearest neighbor plaquette. 

For a 1/6-doped CuO plane we adopt the hypoth- 
esis of charge segregation |2(| into alternating 1/3- 
doped and hole- free Cu-O-... stripes of Goodenough 
et al. 0,0- Under the charge segregation hypothe- 
sis, our findings seem to support the vibronic mecha- 
nism and the traveling charge-density/spin-density wave 
(CDW/SDW) model of Goodenough for the high-T c su- 
perconductivity in cuprates. In this scenario the super- 
conducting state would imply long range "correlation" 
(or phase coherence) among the type of O-atom displace- 
ments shown in Fig. 3(b). We arrive in the end to the 
kind of picture proposed by Egami et al. |21| : electron 
correlation — according to our calculations, the near- 
degeneracy between the quasi-localized ZR like state on a 
CUO4 plaquette and a state whose leading configuration 
is of the form Cud 9 - Op 5 - Cud 9 plus the longer range 
hole hole Coulomb repulsion — is the engine, while the 
driver is the lattice vibration. 

Work combining results of first principles MC SCF cal- 
culations on small clusters and results of calculations for 
a Hubbard type model Hamiltonian on extended CUO4- 
CUO4-... chains is in progress. It is hoped that such 
investigations will enable reasonable estimates of the 
longer range interactions and other (quantitative) predic- 
tions. In a recent study combining the two approaches 
we were able to give a rather realistic description of the 
phase transition in the mixed- valence NaV205 compound 
|29| . "Correlated" calculations for a bi-dimensional Cu- 
O structure to investigate the charge segregation scenario 
in the CuO plane at 1/6-doping are unfeasible at this mo- 
ment. 

The JT/CT effects illustrated by our ab initio results 
are much too robust to be taken just as an artifact of 
our finite clusters or/and of other approximations. We 
believe that these results define a few basic requirements 
that a model Hamiltonian approach should meet for real- 
istic predictions in the optimally doped and under-doped 
copper oxide compounds. 
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